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© Nuclear or electron spin magnetic resonance 
method using multiple radial scans through frequen- 
cy-space. The problem of insufficient sample density 
at the higher image frequencies is solved by deriv- 
ing an edge image, determining a correction image 
therefrom by using a priori knowledge about the 
pixel value distribution and using the correction im- 
age after Fourier transformation to frequency-space. 
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The invention pertains to a method for deriving 
a magnetic resonance image of an object placed in 
a stationary and substantially homogeneous mag- 
netic main field, comprising repeatedly exciting 
said object and sampling magnetic resonance sig- 
nals in the presence of a read gradient superim- 
posed on said main field thereby determining a 
first frequency image of said object on spherical or 
polar coordinates by radially scanning through fre- 
quency-space. 

The invention also pertains to an apparatus to 
perform such method. 

A method as described in the preamble is 
known from the US Patent 4,070,41 1 . In said patent 
is disclosed that magnetic resonance signals are 
induced in an object, notably a patient, present in a 
homogeneous magnetic field and are detected 
while exposing the region of the object to be im- 
aged to a linear gradient superimposed on the 
homogeneous magnetic field. A resonance spec- 
trum is obtained for each of a relatively large 
number of angularly displaced orientations of the 
linear gradient. The image is constructed by a 
mathematical process of back projection. This tech- 
nique uses the fact that the Fourier transform of a 
one-dimensional projection of the spin density re- 
presents a one-dimensional cross-section of the 
three-dimensional Fourier transform of the spin 
density function. 

In the article "NMR Imaging Methods Seen as 
Trajectories in the Reciprocal Space" in Bull. Mag- 
netic Resonance, 6, No. 3, Nov. 1984, p. 140-141, 
is disclosed that the values of the magnetic reso- 
nance signal S(t) fundamentally correspond to the 
values of the Fourier transform of the spin density 
function in frequency-space at positions given by 
dk/dt = -/G, k = (k K ,k y ,k z ) determining the position in 
frequency-space (also called k-space), G =- 
(<3 x ,<3 y ,G z ) determining the magnetic gradient, and 
y being the gyromagnetic ratio. So, in a method as 
described in the preamble, the values of a fre- 
quency image of an image to be produced are 
found as values of the detected signal assigned to 
points along a radial line in frequency-space. 

In US 4,070,61 1 mentioned above is explained 
that radially scanning through frequency-space has 
the distinct disadvantage that information is col- 
lected in an inhomgeneous way, i.e. because of the 
radial scanning the information is more dense near 
the origin of frequency-space, representing the low- 
er frequencies, and there may be insufficient in- 
formation of the higher image frequencies. In US 
4,070,61 1 a method is disclosed in which frequen- 
cy-space is scanned along parallel lines. An im- 
proved version of this method is called spin warp 
and disclosed in US 4,503,222. 

However, especially under circumstances 
wherein the magnetic resonance signal has short 



decay times, which in case of nuclear (proton) 
magnetic resonance depends on the materials or 
tissues examined and which is almost always the 
case with electron spin resonance, it may be ad- 

5 vantageous to start sampling of the magnetic reso- 
nance signal immediately after the end of the ex- 
citation. That is not possible with methods like spin 
warp and spin echo because in those methods an 
echo of the signal is used and time is needed to 

to switch phase encoding gradients. So, especially in 
cases of short decay times, radial scanning may be 
very useful. 

It is an object of the invention to provide a 
radial scan method as described in the preamble in 
75 which the drawback of the less dense information 
acquisition of the higher image frequencies (due to 
radial scanning) is mitigated. 

Another object of the invention is to provide 
such method in which artefacts are greatly re- 
20 duced. 

Still another object of the invention is to pro- 
vide such method in which the probability of the 
image based on its statistical properties is in- 
creased. 

25 These and other objects of the invention are 

achieved in a method as described in the preamble 
by following steps: 

a) determining an edge image in image-space, 

b) deriving from said edge image a correction 
so image in image-space with pixel values derived 

from pixel values of said edge image and a 
predetermined probability function of edge im- 
age values, 

C) determining a corrected frequency image in 

35 frequency-space, 

d) Fourier transforming said corrected frequency 
image to image-space. 
It can be shewn, reference is made to the 
article "Ringing artefact reduction by an efficient 

40 likelihood improvement method" in SPIE Vol. 1137, 
Science and Engineering of Medical Imaging 
(1989) p. 84-90, that an edge image, i.e. an image 
which is derived from a normal image and in which 
mainly the edges are shown between otherwise 

45 relatively uniform areas, has a characteristic dis- 
tribution of pixel (picture element) values. Accord- 
ing to the invention this property is used to deter- 
mine a correction image which is used to correct 
the original image only in so far as required by 

50 insufficient information density. How in that process 
has to be transformed between image-space and 
frequency-space will be explained in the descrip- 
tion with reference to the drawings hereinafter of an 
embodiment of the invention. 

55 In one embodiment of the invention, step a) 

involves 

a1) determining a second frequency image on 
cartesian coordinates from said first frequency 
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image, 

a2) determining an edge frequency image from 
said second frequency image by multiplication 
with an edge function, 

a3) Fourier transforming said edge frequency 
image to image-space obtaining said edge im- 
age in image-space. 
By multiplying the pixel values of a frequency 
image with a factor which is larger according as the 
frequency of a pixel is larger the edges in the 
image (in image-space) become emphasised. Ex- 
periments which led to the present invention have 
shown that a real edge image is obtained by mul- 
tiplication with the edge function j k | a with 0.5 < a 
< 2.0. Best results are obtained for certain images 
with a = 1 .1 . 

in another embodiment of the invention, step a) 
involves 

a1) determining a second frequency image on 
cartesian coordinates from said first frequency 
image, 

a2) Fourier transforming said second frequency 
image to image-space, 

a3) determining a convolution of said Fourier 
transformed second frequency image with an 
edge function obtaining said edge image in im- 
age-space. 

As mathematically speaking multiplication in 
the Fourier domain is equivalent to a convolution in 
the image domain, it is also possible to first trans- 
form the image to image-space and then convert 
the image in image-space to an edge image by a 
convolution with a suitable function. Experiments 
which led to the present invention have shown that 
an excellent edge image is obtained by a convolu- 
tion of the image in image-space with the edge 
function (5 Xl0 -5 x ,i) in the x-direction and with the 
corresponding function in the y-direction and if 
applicable in the z-direction. The pixel values of the 
edge image in image-space are then represented 
by vectors. This convolution is equivalent to sub- 
tracting from every pixel its "positive-x neighbour" 
value. 6 x>i is the Kronecker delta, being 1 only when 
x = i. 

In a further embodiment of the invention step 
c) involves 

d) Fourier transforming said correction image in 
image-space obtaining a correction image in fre- 
quency-space, 

c2) determining a corrected edge frequency im- 
age from said edge frequency image and said 
correction image in frequency-space, 
c3) determining said corrected frequency image 
from said corrected edge frequency image. 
If the edge image in image-space is obtained 
from an edge image in frequency-space the latter 
is directly available for correction with the Fourier 
transform of the correction image to frequency- 



space. 

However, in a another embodiment of the in- 
vention step c) involves 

cl) Fourier transforming said correction image in 
5 image-space obtaining a correction image in fre- 
quency-space, 

c2) Fourier transforming said edge image in 
image-space obtaining an edge frequency im- 
age, 

w c3) determining a corrected edge frequency im- 
age from said edge frequency image and said 
correction image in frequency-space, 
c4) determining said corrected frequency image 
from said corrected edge frequency image. 
15 This embodiment is in particular suitable if the 

edge image is obtained in image-space by a con- 
volution with an edge function because then there 
is no edge image available in frequency-space. As 
a convolution with the function (5 Xi0 -5 Xj1 ) is equiv- 
20 alent to a multiplication in frequency-space by 1- 
exp(brkx), k x being in the range of (-i,^), the trans- 
form in frequency-space back from edge image to 
normal (non-edge) image involves in a 2D case 
multiplication with 

25 

(1 -exp(-i-7rk x ), 1 -exp(-i7rky))/{4-2cos(7rk x )-2cos(7rk y )} 

In an embodiment of the invention said correc- 
tion image in image-space is an image constituting 

30 corrections to increase the probability of said edge 
image in image-space 

This means that the correction image contains 
only the corrections to be applied (with a weight 
factor) to the image to be corrected. Another pro- 

35 cedure would be that the correction image is a 
corrected edge image which is to replace (where 
necessary and with a weight factor) the image to 
be corrected. Mathematically both procedures are 
practically equivalent. 

40 In a preferred embodiment of the invention said 

corrected edge frequency image is a weighted 
combination of said edge frequency image and 
said correction image in frequency-space. 

It would not be an optimal method to apply the 

45 calculated corrections to each pixel in frequency- 
space because the centre of frequency-space has 
been sufficiently dense sampled. Therefore the 
corrections are applied with a weight factor which 
is representative for the measure of correctness of 

50 the uncorrected image. For example, the weight 
factor with which the corrections are multiplied, 
may be low, or even zero, for reliable pixel values 
and high for pixels resulting from interpolation with 
insufficient accuracy. 

55 In an embodiment Of the invention the pixel 

values "e"* of said correction image in image- 
space are determined by e' = b.2a 2 .e/(e 2 + a 2 ), "e" 
being a picture value of said edge image obtained 
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in step a), "a" being a half width half maximum 
value of a probability function of "e" and "b" being 
an empirical factor close to 1 . 

Experiments have shown that the pixel value 
distribution of an edge image in image-space is 5 
very close to a Lorentzian probability distribution p- 
(e) = a/{-7r(e 2 + a 2 )}. Assuming that all the pixels are 
uncorrelated the likelihood of the total image is the 
product of all p(e) and it can be shown that the 
largest increase of this likelihood is given by a w 
correction image (constituting the corrections) 
which is proportional to 2a 2 .e/(e 2 + a 2 ) 

An embodiment of the invention especially di- 
rected to making a two-dimensional image of a 
slice of an object comprises the steps of: ?s 

a) applying a slice selection gradient to said 
main field and simultaneously applying a radio- 
frequency excitation pulse to said object thereby 
exciting spin magnetic dipole moments in a 
slice of said object, 20 

b) applying a read gradient to said main field 
with a gradient direction parallel to said slice, 

c) sampling magnetic resonance signals during 
the application of said read gradient 

d) storing samples as values assigned to points 25 
arranged along a radial line in a frequency- 
space defined on cartesian coordinates k x and 

ky, 

e) repeating steps a) - d) a predetermined num- 
ber of times with different directions parallel to oc 
said slice of said read gradient, thereby obtain- 
ing a first electronic image with values arranged 
along radial lines said in frequency-space, 

f) deriving said magnetic resonance image by 
Fourier transformation, 35 

and is characterised in that between steps e) and 
f): 

aa) all values of said first electronic image are 
processed and multiplied with substantially | k j 
*, with 0,5 < a < 2.0, thereby obtaining an 40 
electronic image E(k x ,k y ) with values arranged 
along lines parallel to k x and k y in said fre- 
quency-space, 

bb) a Fourier transform is applied to the elec- 
tronic image E(k x ,k y ) to derive an electronic in- 45 
termed iate image e(x,y) defined in an image- 
space on cartesian coordinates x and y, 
cc) an electronic correction image e'(x,y) is 
formed from e(x,y) according to the formula 

50 

e'(x,y) = b.2a 2 .e(x,y)/{(e(x,y)) 2 + a 2 } 
wherein "a" is a half width half maximum value of a 
probability function of e(x,y) and "b" is an empiri- 
cal factor close to 1 , 

dd) e'(x,y) is transformed to said frequency- 55 
space by a Fourier transform yielding an elec- 
tronic image E'(k XJ k y ) 



ee) an electronic image 

{E(k x ,k y ) - J fi(k x ,k y ).E'(k x ,k y )}. j k 
wherein /3(k x ,k y ) is a weight factor, is determined as 
an input Image for step f). 

The invention is also directed to an apparatus 
comprising means to perform one of the methods 
described above and in particular to an apparatus 
comprising means to determine an edge image 
from an initial image, means to derive a correction 
image from said edge image, means for correcting 
in frequency-space said initial image with said cor- 
rection image and means for deriving said mag- 
netic resonance image from said corrected initial 
image. 

These and other aspects of the invention will 
be further explained with reference to the accom- 
panying drawings. 
The drawings show in: 

Fig. 1: diagrammatically a magnetic imaging 
apparatus suitable for the method of 
the invention; 
Fig. 2: a sequence of RF and gradient 
pulses suitable for the method of the 
invention; 

Fig. 3: a representation of radial scan in k- 
space; 

Fig. 4: a representation of radial scan in k- 
space; 

Fig. 5: a radial scan image showing some 
artefacts caused by radial scan with 
insufficient resolution of the higher 
image frequencies 

Where in the preceding or following description 
or in the claims the expression "image" is used, 
not only a visual image is meant such as it may be 
observed on a monitor, but also an electronic im- 
age in digital form, being processed or stored in 
memory means. A digital image consists of picture 
elements also called pixels. 

An electronic image may be in image-space 
defined on coordinates x, y and z or in frequency- 
space (also called k-space) defined on coordinates 
k x , k y , and k z . An image in image-space may be 
transformed into an image in frequency-space (and 
vice versa) in a manner well-known to a person 
skilled in the art of magnetic resonance imaging 
(MRI) using a so-called Fourier transform. 

As is customary in this art, functions in fre- 
quency-space will be indicated with capitals and 
functions in image-space will be indicated in small 
print, unless the context requires otherwise. The 
Fourier transformation of a function f in image- 
space to a function F in frequency-space will be 
indicated with 
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and the Fourier transformation of a function F in 
frequency-space to a function f in image-space will 
be indicated with 

A convolution of two functions in image-space or in 
frequency-space will be indicated with <s>, so: f®g 
and F®G respectively. 

Underscored characters generally represent 
vectors. 

As is also customary in the art of magnetic 
resonance imaging, the signals are sampled and 
digitalized and the further image processing is digi- 
tal. All that is implicit in this description. 

The radial scan used in the present invention 
may be used to image an object in three dimen- 
sions (3D) or may be used to image a slice of an 
object in two dimensions (2D). Accordingly refer- 
ence may be made to spherical coordinates (3D) or 
to polar coordinates (2D). The embodiment de- 
scribed hereinafter is a 2D imaging method for the 
ease of the description, but a person skilled in the 
art will have no dificulty in understanding how a 3D 
method works because the additional dimension 
only means more computational effort. 

In the description of an embodiment hereinafter 
reference is made to magnetic resonance, this im- 
plies nuclear resonance which in a specific situ- 
ation is proton resonance and it also implies elec- 
tron spin resonance. 

Magnetic resonance imaging is very useful in 
medical applications, that needs no further explana- 
tion. The object to be imaged is then a patient or 
generally a human being to be examined. 

In Figure 1 a magnetic resonance apparatus 1 
is diagrammatically shown. The apparatus com- 
prises a set of main magnetic coils 2 for generating 
a stationary homogeneous main magnetic field and 
several sets of gradient coils 3, 4 and 5 for super- 
imposing additional magnetic fields with control- 
lable strength and having a gradient in a selected 
direction. Conventionally, the direction of the main 
magnetic field is labelled the z direction, the two 
directions perpendicular thereto the x and y direc- 
tions. The gradient coils are energised via a power 
supply 11. The apparatus further comprises excita- 
tion means 6 for emitting radio-frequency pulses 
(RF-puIses) to an object or body 7, the excitation 
means 6 being coupled to modulating means 8 for 
generating and modulating Of the RF-pulses. Also 
provided are means for receiving the NMR-signals, 
these means can be identical to the excitation 
means 6 or be separate. If the excitation means 
and receiving means are identical, as shown in 
figure 1, a send-receive switch 9 is arranged to 
separate the received signals from the excitation 



pulses. The received magnetic resonance signals 
are input to receiving and demodulating means 10. 
The excitation means 6 and modulating means 8 
and the power supply 11 for the gradient coils 3, 4 

5 and 5 are steered by a control system 12 to 
generate a predetermined sequence of RF-pulses 
and gradient field pulses. The demodulation means 
is coupled to a data processing unit 14, for exam- 
ple a computer, for transformation of the received 

w signals into an image that can be made visible, for 
example on a visual display unit 15. 

If the magnetic resonance apparatus 1 is put 
into operation with an object or body 7 placed in 
the magnetic field, a small excess of magnetic 

15 dipoie moments (nuclear or electron spins) in the 
body will be aligned in the direction of the mag- 
netic field. In equilibrium, this causes a net mag- 
netisation Mo in the material of the body 7, directed 
in parallel with the magnetic field. In the apparatus 

20 1 the macroscopic magnetisation Mo is manipu- 
lated by radiating to the body RE-pulses having a 
frequency equal to the Larmor frequency of the 
dipoie moments, thereby bringing the dipole mo- 
ments in an exited state and re-orienting the mag- 

25 netisation M 0 . By applying the proper RF-pulses, a 
rotation of the macroscopic magnetisation is ob- 
tained, the angle of rotation is called the flip-angle. 
The introduction of variations in the magnetic field 
by applying gradient magnetic fields influences the 

30 behaviour of the magnetisation locally. After the 
application of RF-pulses, the changed magnetisa- 
tion will tend to return to a state of thermal equilib- 
rium in the magnetic field, emitting radiation in the 
process. A well chosen sequence of RF-pulses and 

35 gradient field pulses causes this radiation to be 
emitted as (decaying) magnetic resonance signals 
which provide information about the density and/or 
relaxation times of a certain type of nuclei, for 
example hydrogen nuclei, and the substance in 

40 which they occur. By analysis of the emitted sig- 
nals and presentation of it in the form of images, 
information about the internal structure of the ob- 
ject or body 7 is accessible. The decaying mag- 
netic resonance signal directly after the excitation 

45 is called FID (Free Induction Decay). By proper 
application of the gradients and/or application of 
additional excitation (RF) pulses also an echo of 
the FID may be produced and used to derive 
image information. 

50 The fundamental expression governing mag- 

netic resonance imaging is co = 7 B, meaning that 
the (local) Larmor frequency o> is the gyromagnetic 
ratio y times the (local) magnetic field B. Informa- 
tion about the position is obtained via the gradients 

55 which may be superimposed on the homogeneous 
main field. As mentioned herein before it has been 
shown that the values of the magnetic resonance 
signal S(t) correspond to the values of the Fourier 
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transform of the spin density function in frequency- 
space at positions given by dk/dt = 7G, with k =- 
(kx,ky,k z ) determining the position in frequency- 
space and G = (G x ,G y ,G z ) determining the mag- 
netic gradient. 

For a more detailed description of magnetic 
resonance imaging (MRI) and MRI-devices refer- 
ence is made to the extensive literature on this 
subject, for example to the book "Practical NMR 
Imaging", edited by M.A. Foster and J. M.S. Hutch- 
inson, 1987, IRL Press. 

In Figure 2 a sequence of RF and gradient 
pulses is shown which is used in a 2D magnetic 
resonance method with radial scan. On time axes 
20, 21, 22, 23 and 24 respectively are shown the 
envelope of the RF pulse, slice selection gradient 
G2, gradients G x and G y which together determine 
the amplitude and direction of a read gradient and 
magnetic resonance signal S. One radial scan com- 
prises three time intervals indicated as 1, 2 and 3. 
After a waiting interval 4 the following radial scan 
starts again with an interval 1. The waiting interval 
4 may be left out. 

During interval 1 a slice in the object to be 
examined is excited. The position and thickness of 
the slice is determined by the frequency spectrum 
of the RF pulse 25 of which only the envelope is 
shown and the amplitude of slice selection gradient 
G z . The slice is perpendicular to G z . The known 
rephasing lobe 27 of the slice selection gradient 
(interval 2) serves to achieve a well defined phase 
of all the spins throughout the thickness of the 
Slice. The read gradient formed by gradients G x 
and G y (28 and 29) determines the direction of the 
radial scan, i.e. the direction of k determined by 
angle <#> in Figure 3. During interval 3 the magnetic 
resonance signal 30 is detected and sampled. The 
values of the samples are assigned to positions in 
k-space corresponding to the expression k = 7Gt 
Generally one will take j G j = V(G X 2 + G y 2 )- 
= constant so that each radial scan is passed 
through with the same speed. <f> is determined by 

tgW) =G y /G x . 

Figure 3 is a representation of k-space. A num- 
ber of radial scans 31 are indicated. On radial scan 
32 schematically the positions of some measured 
samples are shown. Each radial scan may start 
with a very short scan in the opposite direction in 
order to better detect the very beginning of the free 
induction decay in the form of an echo. Although 
the radial scans are shown equidistant (the same 
A(f> between the scans) this is not essential. If the 
object under investigation is oblong then it may be 
advantageous to have a higher sample density in 
the length direction, i.e. a smaller A<j> in the length 
direction of the object then in the direction perpen- 
diculat thereto. 



Figure 4 is again a representation of k-space. 
As will be apparent, there are areas 41 of larger j k 
j =V(k x 2 + ky 2 ) where k-space is insufficiently filled 
with data. As large [ k j corresponds to the higher 

5 frequencies in the image, the reproduction of the 
finer image details is in jeopardy. Supposing the 
unknown data are zero causes artefacts in the 
image displayed. 

After k-space has been completely scanned, 

70 the data as a function of k determined by coordi- 
nates | k j and £ are transformed to cartesian 
(rectangular) coordinates k x and k y . A simple ap- 
proach would be interpolation but more sophisti- 
cated methods are available and give much better 

15 results. In the embodiment here described a so- 
called gridding algorithm is used which implies 
convolution with a sine function followed by resam- 
pling. One such gridding algorithm Is described in 
the article "A Fast Sine Function Gridding Algo- 

20 rithm for Fourier Inversion in Computer Tomog- 
raphy" in IEEE Transactions on Medical Imaging, 
Vol. MI-4, No. 4, December 1985, p.200-207. 

After gridding the values now on k x , k y coordi- 
nates are multiplied with j k | " with a -1.1. The 

25 value of a has to be determined empirically and 
may lie between 0.5 and 2.0. Also in a 3D method 
a need not to be much larger than 1 although 
intuitively one would expect a higher value. The 
purpose is to obtain an image which after Fourier 

30 tranformation to image-space is an edge image, i.e. 
an image in which only the edges between rela- 
tively uniform areas are shown. As determining an 
edge image is a kind of differentiation of the image 
and in one dimension differentiation is equivalent to 

35 multiplication with k in k-space, one may expect 
that multiplication with some power of j k j may 
result in a usable edge image in image-space. 
Gridding and multiplication with [ k | not nec- 
essarily but preferably in that order, yield an edge 

40 image E(k x ,k y ) in k-space and after transformation 
to image-space an edge image e(x 5 y). 

It is known that such an edge image like e(x,y) 
has a characteristic distribution of pixel values. This 
distribution is not Gaussian as one would expect 

45 perhaps in first instance but more Lorentzian. See 
e.g. the article about ringing artefact reduction in 
SPIE Vol. 1137 (1989) mentioned herein before. 
However, because the edge image obtained is 
based on insufficient data (caused by the radial 

50 scan) it contains all kinds of artefacts. On the other 
hand, the amplitude of these artefacts appears to 
be relatively small because the center of k-space 
was sufficiently sampled. As, because of the mul- 
tiplication with | k | a in k-space, the uniform areas 

55 of e(x,y) have relatively low pixel amplitudes and 
the edges have relatively high pixel amplitudes, 
there is a range of pixel amplitudes which are less 
likely to belong to a correct edge image and more 
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likely to belong to undesired artefacts. Experiments 
have shown that suppressing to a certain extent 
intermediate pixel values is an effective method to 
correct an edge image affected with such artefacts. 
Analogous to said ringing artefact reduction a cor- 
rection image e'(x,y) is determined from e(x,y) ac- 
cording to the formula 

e'(x,y) = b.2a 2 .e(x,y)/{(e(x,y)) 2 + a 2 } 

In this expression is "a" the half width half 
maximum value of the Lorentzian probability dis- 
tribution and "b" an empirical factor close to 1 . "a" 
may depend on the type of image and the way the 
edge image is derived from it. Good results have 
been obtained with "a" close to 1.7 time the aver- 
age value of the absolute values of e(x,y). 

However, it would not be an optimal method if 
the corrected edge image would be simply ob- 
tained as e(x,y)-e'(x,y). According to the invention 
use is made of the understanding that the Fourier 
transform of e'(x.y) to k-space, E^k^ky), ideally 
should be zero on positions in k-space with suffi- 
cient sampling density. This ideal is not reached, 
therefore, e'(x,y) is first transformed to k-space 
yielding E^k^ky). Then, in k-space, E*(k x ,k y ) is mul- 
tiplied with a weight factor /8(k x ,k y ) and the product 
is subtracted as a correction image from E(k x ,k y ). 
This result is multiplied with j k | ~" in order to go 
back from edge image to normal image in k-space. 
This final image in k-space is then Fourier trans- 
formed to image-space in order to obtain the de- 
sired improved image. 

The weight factor /S(k XJ k y ) can be determined in 
various ways. One can take j9-0 in areas of k- 
space with sufficient sampling density and = 1 
elsewhere. Known gridding methods also provide 
information about the density of neighbouring mea- 
sured data. Then one can determine an auxiliary 
array on k x ,k y coinciding with the image data points 
and containing a reliability factor w(k Xl k y ) which 
approaches 1 in areas with sufficient sampling den- 
sity and approaches 0 in areas where all sampling 
points are far away. Good results have been ob- 
tained by taking 

P(KM y ) =0.5 4 /(w(k X( k y )) 4 +0.5 4 ) 

as a continuous weight factor. 

Artefacts due to "missing data" are not nec- 
essarily removed in a single iteration. One of the 
reasons is that the artefacts in the uncorrected 
image may be so large that they are seen, with 
relatively high likelihood, as belonging to actual 
image structures. In that case, the first iteration 
only manages to slightly reduce the artefacts. Ex- 
periments leading to the present invention have 
shown that this slight improvement of the distinc- 



tion between artefacts and real edges makes it 
possible to improve the correction with the applica- 
tion of (a) following iteration(s). 

Another promising way to use the present in- 
5 vention Is to present the edge values e(x,y) by a 
vector instead of scalar value. This vector e(x,y) = - 
(e x ( x >y)>ey( x .y)) represents the image intensity gra- 
dient on every location in image-space. The Fourier 
transform of such a vector 

10 

^{e(x,y)}=E(k„k y ) 

can be seen as 

15 

E(k > ,k y )=C^{e x (x,y)},^ fe y (x,y)}). 

This embodiment of the invention can be seen as 
20 replacing the multiplication with j k | * by a mul- 
tiplication with k = (k x ,k y ). The inversion of the mul- 
tiplication with k is a multiplication with k/ [ k j 2 . 
The correction image is then determined by 

25 e'(x,y) = b.2a 2 .e(x,y)/{ | e(x,y) | 2 + a 2 } . 

These and other modifications of the methods 
described make use of the essence of the inven- 
tion. 

30 Figure 5 shows a 2D transaxial head scan 

obtained with radially scanning through frequency- 
space along 180 radii (2° intervals). Apparent are 
socalled streaks radiating from structures in the 
image. These artefacts can be removed to a large 

35 extent if the present invention is used. Because the 
artefacts are hardly detectable on a low resolution 
medium as a photocopy a corresponding image 
with reduced artefacts is not shown. 

Finally, it should be remarked that in an elec- 

40 tron spin resonance (ESR) imaging method with 
relaxation times of the order of 1 us slice selection 
techniques are not appropriate. Then a 3D radial 
scan method often will be used and a method 
according to the invention will be very useful. 

45 

Claims 

1. Method for deriving a magnetic resonance im- 
age of an object placed in a stationary and 

so substantially homogeneous magnetic main 

field, comprising repeatedly exciting said ob- 
ject and sampling magnetic resonance signals 
in the presence of a read gradient superim- 
posed on said main field thereby determining a 

55 first frequency image of said object on spheri- 

cal or polar coordinates by radially scanning 
through frequency-space, characterised by the 
steps of: 
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a) determining an edge image in image- 
space, 

b) deriving from said edge image a correc- 
tion image in image-space with pixel values 
derived from pixei values of said edge im- 5 
age and a predetermined probability func- 
tion of edge image values, 

c) determining a corrected frequency image 
in frequency-space, 

d) Fourier transforming said corrected fre- w 
quency image to image-space. 

2. Method as claimed in claim 1, characterised in 
that step a) involves: 

a1) determining a second frequency image 1& 
on cartesian coordinates from said first fre- 
quency image, 

a2) determining an edge frequency image 
from said second frequency image by mul- 
tiplication with an edge function, 20 
a3) Fourier transforming said edge frequen- 
cy image to image-space obtaining said 
edge image in image-space. 

3. Method as claimed in claim 1, characterised in 25 
that step a) involves: 

a1) determining a second frequency image 
on cartesian coordinates from said first fre- 
quency image, 

a2) Fourier transforming said second fre- 30 
quency image to image-space, 
a3) determining a convolution of said 
Fourier transformed second frequency im- 
age with an edge function obtaining said 
edge image in image-space. 36 

4. Method as claimed in claim 2, characterised in 
that step c) involves: 

c1) Fourier transforming said correction im- 
age in image-space obtaining a correction 40 
image in frequency-space, 
c2) determining a corrected edge frequency 
image from said edge frequency image and 
said correction image in frequency-space, 
c3) determining said corrected frequency 45 
image from said corrected edge frequency 
image. 

5. Method as claimed in claim 3, characterised in 

that step c) involves: 50 
c1) Fourier transforming said correction im- 
age in image-space obtaining a correction 
image in frequency-space, 
c2) Fourier transforming said edge image in 
image-space obtaining an edge frequency 55 
image, 

c3) determining a corrected edge frequency 
image from said edge frequency image and 



said correction image in frequency-space, 
c4) determining said corrected frequency 
image from said corrected edge frequency 
image. 

6. Method as claimed in any one of claims 1 - 5, 
characterised in that said correction image in 
image-space is an image constituting correc- 
tions to increase the probability of said edge 
image in image-space. 

7. Method as claimed in claim 4, 5 or 6, charac- 
terised in that said corrected edge frequency 
image is a weighted combination of said edge 
frequency image and said correction image in 
frequency-space. 

8. Method as claimed in claim 1 , characterised in 
that pixel values "e'" of said correction image 
in image-space are determined by e' = b.2a 2 .e/- 
(e 2 + a 2 ), "e" being a picture value of said edge 
image obtained in step a), "a" being a half 
width half maximum value of a probability 
function of "e" and "b" being an empirical 
factor close to 1 . 

9. Method for deriving a magnetic resonance im- 
age of an object placed in a stationary and 
substantially homogeneous magnetic main 
field, comprising the steps of: 

a) applying a slice selection gradient to said 
main field and simultaneously applying a 
radio-frequency excitation pulse to said ob- 
ject thereby exciting spin magnetic dipole 
moments in a slice of said object, 

b) applying a read gradient to said main 
field with a gradient direction parallel to said 
slice, 

c) sampling magnetic resonance signals 
during the application of said read gradient, 

d) storing samples as values assigned to 
points arranged along a radial line in a 
frequency-space defined on cartesian co- 
ordinates k x and k y , 

e) repeating steps a) - d) a predetermined 
number of times with different directions 
parallel to said slice of said read gradient, 
thereby obtaining a first electronic image 
with values arranged along radial lines said 
in frequency-space, 

f) deriving said magnetic resonance image 
by Fourier transformation, 

characterised in that between steps e) and f): 
aa) all values of said first electronic image 
are processed and multiplied with substan- 
tially j k |* with 0.5 < ot < 2.0, thereby 
obtaining an electronic image E(k XJ k y ) with 
values arranged along lines parallel to k x 
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and k y in said frequency-space, 
bb) a Fourier transform is applied to the 
electronic image E(k x ,k y ) to derive an elec- 
tronic image e(x,y) defined in image-space 
on cartesian coordinates x and y, 5 
cc) an electronic correction image e'(x.y) is 
formed from e(x,y) according to the formula 

e'(x,y) = b.2a 2 .e(x,y)/{(e(x,y)) 2 + a 2 } 
wherein "a" is a half width half maximum value io 
of a probability function of e(x,y) and "b" is an 
empirical factor close to 1, 

dd) e'(x,y) is transformed to said frequency- 
space by a Fourier transform yielding an 
electronic image E^k^ky) 15 
ee) an electronic image 

{E(k x ,k y ) - iS(k x ,k y ).E'(k x .ky)}. | k I 
wherein /3(k x ,k y ) is a weight factor, is deter- 
mined as an input image for step f). 20 

10. Apparatus for deriving a magnetic resonance 
image, characterised by means to perform a 
method as claimed in claims 1 - 9. 

25 

11. Apparatus as claimed in claim 10, characteris- 
ed by means to determine an edge image 
from an initial image, means to derive a cor- 
rection image from said edge image, means 

for correcting in frequency-space said initial 30 
image with said correction image and means 
for deriving said magnetic resonance image 
from said corrected initial image. 
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